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ABSTRACT 

In this paper we present a procedure able to identify and reject showers with the 
core outside the ARGO-YBJ carpet boundaries. The efficiency of this method 
is investigated for different primary energies and fiducial areas. A comparison 
of the results for gamma and proton induced showers is also presented. 



1 Introduction 

The ARGO-YBJ detector, currently under construction at the Yangbajing Lab- 
oratory (P.R. China, 4300 m a.s.l.), is a full coverage array of dimensions 
~ 74 X 78 realized with a single layer of RPCs. The area surrounding 
this central detector (carpet), up to ~ 100 x 110 m^, is partially 50%) in- 
strumented with other RPCs. The basic element is the logical pad (56 x 62 cm^) 
which defines the time and space granularity of the detector. The layout of the 



detector is shown in Fig. ^ The detector is subdivided in 6 x 2 RFC units 
(Clusters, the rectangles of Fig. P). 

Showers of sufficiently large size will trigger the detector even if their 
core is located outside its boundaries. The corresponding core positions are 
generally reconstructed not only near the carpet edges but also well inside 
them. As a consequence, sofisticated algorithms to reduce the contamination 
of external events are needed. The goal is to identify and reject a large fraction 
of external events before exploiting any reconstruction algorithm only by using 
some suitable parameters. 

The rejection of external events is important because a large difference 
between the true and the reconstructed shower core position may lead to a 
systematic miscalculation of some shower characteristics, such as the shower 
size. Moreover, an accurate determination of the shower core position for se- 
lected internal events is important to reconstruct the primary direction using 
conical fits to the shower front, improving the detector angular resolution or to 
performe an efficient gamma/hadron discrimination. 

In this paper we present a reconstruction procedure able to identify and 
reject a large fraction of showers with core outside the ARGO-YBJ detector. 
The efficiency of this procedure is investigated both for gamma and proton 
induced showers. 

2 Identification of external events 

To perform these calculations we have simulated, via the Corsika code , 
7-induced showers with a Crab-like spectrum E^'^-^) and proton events 
with ~ E^"^-^^, both ranging from 100 GeV to 50 TeV. The 7-rays have 
been simulated for different zenith angles (< 40°), following the daily path 
of the source in the sky. The detector response has been simulated via a 
GEANT3-based code. The core positions have been randomly sampled in an 
area, energy-dependent, large up to 800 x 800 w? centered on the detector. 

In Fig. 121 we show the fraction of showers with the core truly inter- 
nal (external) to a fiducial area approximately equal to the carpet dimensions 
{A fid = 80 X 80 m?) as a function of the number of pads fired on the ARGO- 
YBJ carpet. The upper plot refers to 7-induced showers, the lower to proton- 
induced events. As expected, for low pad multiplicities the events are mainly 
external: about 40% of the 7-induced showers with a multiplicity of 100 - 150 
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Figure 1: The ARGO-YBJ detector. Figure 2: Fraction of events with the 

core IN (OUT) Af.^ = 80 x 80 m^. 



fired pads are external to A fid- But also for higher multiplicities the percentage 
of external events is consistent (about 30%). 

Various parameters based on particle density or time information are un- 
der investigation to identify showers with core position outside a given fiducial 
area 13. In this paper we discuss the performance of the following ones: 

• Position of the cluster with the highest particle density. 

• Position of the cluster row/column with the highest particle density. 

• Mean distance Rp of all fired pads to the reconstructed shower core po- 
sition. 

As an example, in Fig. Owe show the distributions of the positions of the 
cluster with the highest particle density for 7-induced showers which fire more 
than 100 pads on the central carpet. In the plots we compare the events with 
the core truly external to a 80 x 80 ip? fiducial area (solid histograms) and 
the truly internal ones (dashed histograms). To investigate the discrimination 
power of this particular parameter we have simulated a detector completely 
instrumented up to ~ 100 x 110 m^, i.e., containing 14 x 17 clusters. There- 
fore, the cluster coordinates nm from 1 to 14 (X view) and from 1 to 17 (Y 
view) starting from the lower left corner of the carpet (see Fig. The clear 



difference between the IN and OUT showers suggests to tag as external the 
events with the highest particle density in the outer clusters and to reject them 
before exploiting any reconstruction algorithm. 

The mean lateral spread of the shower can be expressed as: 

P T,Zi \rc-n\n, 

where N is the number of fired pads, rc is the reconstructed core position, 
is the position of the i-th fired pad, rii is the number of detected electrons in 
the i-th pad. The Rp distribution for showers reconstructed inside a 80 x 80 
m^ fiducial area is shown in Fig. 0] (solid histogram). The dashed line refers to 
truly IN events while the dotted histogram refers to OUT showers erroneously 
reconstructed as internal. The shower cores have been calculated by means 
of the simple center of gravity method. As can be seen, the parameter Rp 
identifies quite well the events with core outside the carpet. Large distances 
between the true and the reconstructed shower axis lead to larger Rp values. 
This fact offers the possibility to define a cut in Rp to identify these events. A 
conservative choice is to reject showers with Rp > 25 m. 

From these studies it follows that the identification of a large fraction 
of external events can be achieved by defining a suitable fiducial area and a 
combination of cuts in the parameters discussed above. 

3 Maximum Likelihood Method (LLF) 

Different algorithms have been investigated to reconstruct the shower core po- 
sition in the ARGO-YB J experiment ^ . The most performant one is the 
Maximum Likelihood Method. 

If < mi > is the average particle number expected on the i-th pad, then 
the probability of finding Ni particles is 

Therefore Pi{0) — e~<™'> is the probability of finding particles, and Pi{> 
0) = I — e^<™*> is the probability of finding 1 or more particles. The Likelihood 
Function (LF) is given by: LF ~ HiPi- The natural Logarithm of LF (LLF) 
becomes a sum: 



LLF = ln{n,P,) = T,ilnP^{0) + T,jlnPj{> 0) 



(3) 
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Figure 3: Coordinate distributions of 
the cluster with the highest particle 
density for "/-induced events with pad 
multiplicity Nhu > 100. 
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Figure 4: Distributions of the pa- 
rameter Rp for ^-induced events with 
NMt > 100. 



where the index i runs on the not fired pads, while the index j refers to the 
fired pads. Exploiting the relation < mi >= Spad ■ Pii^^^ obtain 

- LLF = Spad ■ ^iPi - T.jln{pj) - Npad{> 0) • ln{Spad) (4) 

where pi — f{Ri/R'i^) ■ -^if is the particle density expected on the i-th pad 
at a distance Ri from the core. The lateral structure function f{Ri/R'j^j) has 
been calculated for the ARGO-YBJ detector: 

with the normalization factor defined by 

C ^(^) (6) 

2ti{R\j)^--^T{A)T{B ~ A) ^ ^ 

where — i?Af/3.944, _Rj\/ being the Moliere radius (133 m at YBJ altitude), 
and A = 1.826, B = 2.924, C = 0.613 . Spad is the pad area and Npad{> 0) is 
the total number of pads fired by the shower. The minimum value of -LLF is 
then chosen as the best fit for the freely varying parameters {xc, yc, A^e}, being 



{xcUc} the core coordinates and the shower size. In the following we will 
refer to this approach as to the "LLFl method". 

We point out that expression Q for -LLF refers to the case of a Poisson 
distribution in which the pads are not fired with probability Pi{0) or fired with 
probability Pi{> 0) = 1 — Pi{0). In our study almost always the fired pads 
have particle multiplicity 1, and therefore such a simple discrimination can be 
made. However, if we consider a larger area as the whole RPC, the multiplicity 
can also be > 1, and the proper Poisson distribution on the fired RPCs appears 
more adequate. In this case the sum on fired elements is: 



-j:jlnPj{> 0) = --J:jNjln{pj)~ln{SRPc)^jNj+^jln{Nj\) + SRPc^jPj (7) 



where pj is the particle density expected on the j-th RPC at a distance Rj from 
the core, Nj is the recorded particle number and Shpc is the RPC area. The 
shower size can be calculated via the equation 



We call this calculation the " LLF2 method" . As a consequence, we expect 
that the differences between LLFl and LLF2 increase, for a fixed area, with the 
particle density. In Fig. we compare the shower core position resolution 



Ccore — y'^x + Calculated by applying the LLFl and LLF2 methods on the 
RPCs for 7-induced showers with the core randomly sampled inside a 80 x 80 
area. As expected, the resolution worsens with the multiplicity if the LLFl 
approach is applied when the number of particles hitting the RPC is > 1. 
We note that for very low multiplicities {Nhu < 80) method LLFl is more 
performant than LLF2. Indeed, the algorithm based on the RPC occupancy 
(LLFl) provides a better representation of the hit distribution in very poor 
showers. For very high multiplicities {Nhn > 10^) the shower core position is 
determined by LLF2 with an uncertainty < 1 m. 

In Fig. 1^ cTcore obtained with the LLF2 method is shown as a function 
of multiplicity for 7 and proton induced showers with core randomly sampled 
in an area, energy-dependent, large up to 800 x 800 m^. The fiducial area 
is 80 X 80 171^. The increase of the sampling area worsens the performance 
of the method because of the contamination by external events (see Fig. 
This effect can be limited exploiting the cuts described in section [3 As an 
example, in Fig. \7\we compare the shower core position resolutions for 7 and 
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Figure 5: Comparison between the 
shower core position resolutions ob- 
tained with LLFl and LLF2 methods. 
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Figure 6: Shower core position reso- 
lutions obtained with LLF2 for 7 ( cir- 
cles) and proton (triangles) showers. 



proton showers after the following selection procedure: (1) Rejection of the 
events whose highest density clusters are on the guard ring {X = {1, 14}; Y = 
{1, 17}) or on the boundaries of the central carpet {X = {3, 12}; Y — {3, 15}); 
(2) Rejection of the events whose highest total density rows or columns are 
respectively in positions {1, 3, 15, 17} or {1, 3, 12, 14}; (3) Reconstruction of core 
coordinates {Xc, Yc} using LLF2; (4) Further rejection of events with Rp > 25 
m. 

For events with a multiplicity of 100 — 150 hits (Jcore improves from 6.9 
to 4.4 m for 7-induced showers and from 24.6 to 21.3 m for proton events. For 
any given multiplicity the shower core position resolution is better for 7 than 
for proton induced showers, due to the larger lateral particles spread in the 
latter. 

4 Results 

The procedure (1) - (4) is one of the possible procedures to reject external 
events in the ARGO-YBJ detector. In Fig. |Slthe fraction of rejected events 
(internal and external to A fid = 80 x80 , respectively) is reported for 7 and 
proton-induced showers. As can be seen, this algorithm is able to identify and 
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Figure 7: Comparison between the 
shower core position resolutions for 7 
(circles) and proton (triangles) events 
after the selection procedure (1) - (4)- 
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Figure 8: Fraction of truly internal 
(Fin) and external (Font) events re- 
jected by the selection procedure (1) - 

U). 



reject a large fraction of external events: for photons this percentage is always 
larger than 95%. For low multiplicities {Nhu < 150) a significative fraction of 
internal events is erroneously rejected, especially in proton-induced showers. 

The fraction of internal and external events rejected with the above pro- 
cedure is shown for different fiducial areas, separately for 7-induced showers 
(Fig. O and for proton-induced events (Fig. I1U|I . The fraction of discarded 
showers increases with Afi^ mainly at high multiplicities. For 7-showers more 
than 94% of external events are rejected for the investigated fiducial areas. The 
fraction of internal events erroneously rejected is less than 30% for high mul- 
tiplicities {Nhit > 150). The performance of the procedure is lower for proton 
showers. 

Due to the larger lateral particles spread in proton showers the 7/hadron 
relative trigger efficiency is smaller for external events than for internal ones. 
Thus, a suitable choice of the fiducial area A fid may improve the overall detec- 
tor sensitivity to point 7-ray sources even if the total effective area is decreased. 
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Figure 9: Fraction of truly internal 
(Fin) and external (Font) 7 showers 
rejected by the selection procedure (1) 
- (4) for different fiducial areas. 



Figure 10: Same as Fig. but for 
proton showers. 



5 Conclusions 

The fraction of external events which triggers the ARGO-YBJ detector is rel- 
evant, mainly for low multiplicities (see Fig. EJ. In this paper we discussed 
a possible procedure to tag and reject a relevant fraction of external showers 
before the calculation of the shower core position via the LLF2 method, thus 
saving CPU time. We applied this algorithm to a sample of 7-induced showers 
with a Crab-like energy spectrum. For comparison, we studied the perfor- 
mance of the procedure for proton-induced showers with an energy spectrum 

With a suitable choice of the fiducial area and applying simple cuts we 
can reject more than 94% of the external 7 showers and more than 80% of 
the external proton events, saving more than 50% of the internal events (for 
Nhit > 70) . As a consequence, we expect an improvement in the reconstruction 
of some shower characteristics, such as the primary direction 
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